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Following the intriguing hints of deviations from the Standard Model 
in rare B meson decays, searches for rare and forbidden decays of charm 
hadrons become a hot topic again. We present recent results on Flavour 
Changing Neutral Current —)■ pCpT, p'^p' -TT+TT , Dj) ^ 

Tr^p~^p~ and LFV —)■ n^p^p^ obtained at LHCb. 

Some future prospects on the field of charm rare decays are also ex¬ 
plained as well as some predictions on what is expected after LHC Run 
H and LHCb Upgrade. 
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1 Introduction 


Charm rare decays represent a unique chance to investigate Flavour Changing Neutral 
Currents (FCNC) processes mediated by up-type quarks. These studies, which are 
complemetary to those on B and K sectors, together with those on Lepton Number 
Violating (LNV) processes, could allow to test Standard Model predictions and have 
a privileged point of view on confirming or rejecting New Physics (NP) theories. In 
fact, FCNC processes are highly suppressed on Standard Model and only allowed at 
loop level and the Glashow-Iliopoulos-Maiani (GIM) suppression affects D decays, 
involving c —)■ processes, more than B ones, due to the absence of a high-mass 

down-type quark. 

The extremely low branching fraction predicted by the SM make the D —)■ 
decays be perfect candidates to investigate NP models predicting enhancements up 
to several orders of magnitude. 

Furthermore, multibody semileptonic decays allow angular studies such as searches 
for Forward-backward asymmetries, which could also be enhanced by some NP effects 
to 0(1%) d] [2]. 

2 The LHCb detector 

The Large Hadron Collider-beauty (LHCb) detector is a single-arm forward spec¬ 
trometer, aimed to study b- and c-hadrons rare decays, CP violation, test the quark 
model and investigate the physics beyond the Standard Model (matter-antimatter 
asymmetry). For this aim, it works under low luminosity conditions, with few p-p 
interactions per bunch crossing to ensure a better reconstructibility of events. These 
properties, together with an excellent muon identihcation, a high momentum resolu¬ 
tion (0.4% < dp/p < 0.6%), very good performance in reconstruction of vertices and 
a high performance trigger, flexible and conhgurable make LHCb be a very suitable 
detector for studying charm rare decays. Moreover, bb/cc pairs are predominately 
produced at high \p\ and 5(2)-10^^ have been produced in LHCb acceptance 

(1.9 < p <4.9) in 3fb~^ of integrated luminosity at y/s = 7-8 TeV. 

3 Search for the rare decay ^ 

The measurement on —)■ pT branching fraction |3] has been performed by 

analysing a data sample corresponding to an integrated luminosity of 0.9 fb~^ of pp 
collisions collected at a centre-of-mass energy of 7 TeV by the LHCb experiment. 

This FCNC process is expected to be extremely rare in SM, due to an additional 
helicity suppression. The short distance contributions (0(10“^®)) [1] are negligible 
with respect to long distance ones which are dominated by a two-photons intermediate 
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state. The current upper limit on B{D^ —)■ 77 ) is 2.6 ■ 10“® [5] and translates into an 
upper bound for SM predictions of 0(10“^^). 

The previous limit on this branching fraction was set by the Belle Collaboration 
to 1.4 ■ 10“^ [^. 

The search for the decay is performed using D*~^ —decays, with the 
D*~^ produced directly at a p-p collision primary vertex. 

The D*^ —)■ D^{7r^7T~)7r^ has been selected as a normalisation mode and the 
branching fraction has been obtained from: 

^ p+/i”) = • — ■ f3{D^ TT+TT-). 

Some control channels {D*~^ 7 r+) 7 r+, —)■ K~n^ without the D*+ 

assumption, J/T —have been used to determine muon identihcation and 
trigger efficiencies and in particular the K~tt^ sample has been used to control 

the misidentification rate of the pion as muon, which was one of the major aspects of 
the analysis. 




Figure 1: (a) Invariant mass difference with in the range 1820-1885 MeV/c^ and (b) invariant 
mass TO^+^-, with Am^+^- = in the range 144-147 MeV/c^ for £)*+ —>• iD°(/r+/r“)7r+ 

candidates. The projections of the two-dimensional unbinned extended maximum likelihood fit are 
overlaid. The curves represent the total distribution (solid black), the D*'^ —)■ (dashed 

red), the combinatorial background (dashed yellow), the Z3*+ —>■ D^{K~'k^)-k^ (dash-dotted blue), 
the D*^ —>• Z?°(7r“^+r'^)7r+ (dash-dotted purple) and the signal £)*+ —^ I?°(^+/r“)7r+ (solid green) 
contribution. 

The backgroud originates from two sources: 

• peaking background (2- or 3-body decays in which hadrons are misidentified 
as muons), reduced by applying tight particle identification criteria; 
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• combinatorial background, due to semileptonic decays of beauty and charm 
hadrons, suppressed by a multivariate selection based on a boosted decision 
tree whose input variables are: the pointing angle xjp of vertex 
and the two muons, minimum of the two muons, Xkf (Kalman hlter) of the 
constrained £t, positively-charged muon angle in the rest frame with respect 
to the flight direction and angle in the D*~^ rest frame with respect to 
the D*~^ flight direction. 



Figure 2: Confidence levels (CLs) (solid line) as a function of the assumed —>■ branching 

fraction and median (dashed line), Ict and 2cr bands of the expected CLs, in the background- 
only hypothesis, obtained with the asymptotic CLs method. The horizontal line corresponding to 
CLs=0.05 is also drawn. 

The result obtained for the upper limit on this branching fraction is: 

B{D^ fi+p-) < 6.2(7.6) ■ 10-9 at 90% (95%) CL 

and represents an improvement by a factor 20 with respect to the previous limit, 
but still lets 2 orders of magnitude to be investigated before reaching the value pre¬ 
dicted in the SM. 


4 Search for —)► tt+tt decay 

The —)■ 7r~^F ~decay branching fraction is dominated by long distance con¬ 
tributions such as V{fi~^p~)7r~^F~ (0(10“®)), in which muons come from a 

resonance (p, 0, p). The measurement performed at LHCb |9] aims at constraining 
the short distance contribution to this decay. This is an FCNC process whose branch¬ 
ing fraction is expected to be lower than lO”® in the SM and for this reason it can be 
used to probe NP. 

The best limit before this result was 3.0 ■ 10“® at 90% CL (conhdence level) and 
was performed by E791 Collaboration [7]. 
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The data sample used corresponds to an integrated luminosity of 1.0 fb~^ at y/s 
= 7 TeV collected by the LHCb detector. 

The sample of D^s is originated from the decay D*~^ —Zl°7 r’'‘.The analysis is 
performed in di-muon mass regions, to minimize the leakage coming from resonances. 
The so-called signal regions are dehned as 250 < < 525MeV/c^ and > 

llOOMeV/c^. 

The selection has been performed by combining a multivariate analysis based on 
BDT(input variables:6 'd, of decay vertex and flight distance, yfp, p and pt of 
all hnal state tracks, of the vertex and pt of the D*^ candidate, maximum distance 
of closest approach between all pairs of tracks forming the and D*~^ candidates) 
and some requirements on muon particle identihcation. A peaking background is 
originated from the misidentihcation of pions in the decay —)■ 7r’''7r“7r’''7r“, while a 
secondary source is combinatorial one. 

The measurement has been normalised to a reference sample of (^{pi^pr) 

decays. Unfortunately, the branching fraction of this decay has not been measured di¬ 
rectly. Its value has been derived from the £t fractions of the various 
contributions found by the amplitudes analysis for the —)■ decay 

reported in [8], and from the ’^citio. This procedure results in a system¬ 

atic uncertainty of 17% on the branching fraction of the normalisation mode, which 
dominated the total uncertainty: —)■ 7r+7r“0(/r+/i“)) = (5.2 ± 1.1) ■ 10“’^. 

The —>■ yields have also been fitted in the regions where the 0 

and p resonances peak. In order to evaluate the size of the leakage from resonances 
into the signal regions, the fraction of such decays lying in these regions was esti¬ 
mated assuming Breit-Wigner lineshapes. No interference between amplitutes has 
been considered, since the corresponding effects are negligible at the level of precision 
reachable. 

The signal and background yields have been measured using an unbinned maxi¬ 
mum likelihood £t of two-dimensional [mjr 7 r/i^,Am] distributions in ranges 1810-1920 
MeV/c^ and 140-151.4 MeV/c^ respectively, with Am=m,r 7 r/i^( 7 r — (Fig|^. 

The result obtained for the branching fraction in the signal regions has been 
extrapolated to the whole dimuon mass window by assuming a phase space model. 

The upper limit for the branching fraction is set to: 

B{D^ TT+TT-p+p-) < 5.5(6.7) ■ 10-^ at 90% (95%) CL 

This limit improves the previous one by a factor almost 70, but still lays 2 orders 
of magnitude above SM predictions. 
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Figure 3: [1] Distribution of m(7r+7r“/i+/i“) for —>■ 7r+7r“^+/j,“ candidates in the a) low 

m(^+/i“), h)p/uj, (c) (/), and (d) high-m(^+/r“) regions, with Am in the range 144.4-146.6 MeV/c^. 
[2] Distribution of Am for — )■ 7r+7r“/i+/r“ candidates in the a) low m(/r+/x“), h)p/uj, (c) (f), and 
(d) high-m(/r’''/r“) regions, with invariant mass in the range 1840-1888 MeV/c^. The data are 
shown as points (black) and the t result (dark blue line) is overlaid. The components of the t are 
also shown: the signal (lied area), the —)■ 7r+7r“7r+7r“ background (green dashed line) and the 
non-peaking background (red dashed-dotted line). 


5 Search for —)► 7r+/i+/i / —)► tt decays 

The measurement of the —)■ and D — )■ [ID] branching fractions 

has been performed at LHCb using proton-proton collision data, corresponding to an 
integrated luminosity of 1.0 fb~^ at ^/s=7 TeV recorded in 2011. 

While the —)■ decay is dominated by D~^ —)■ 7 r+l/(/r+/i“) long distance 

contributions and also receives contribution from the FCNC process, we aimed at 
probing that the Df —proceeds via a weak annihilation CKM suppressed. 
The latter can also be used to normalise a potential weak annihilation contribute in 
the D'^ —vr^/i■*■/!“. The limits set for these processes before LHCb measurement 
were 

B{D^ 7r+/x+/i") < 3.9 ■ 10“® [12] 

B{D^ —>■ Tr^/x+zi”) < 2.6 ■ 10“® [T3| 

The —)■ decay is a LNV (Lepton Number Violation) process, and 

thus prohibited in the SM, since it could only happen via lepton mixing mediated by 
non -SM particle as a Majorana neutrino. The previous limits established for these 
decays were 0(10“®) for the D~^ —)■ 7 r“/x+/i+ [TT] and one order of magnitude above 
for the Df —>■ tt~ [TT] . 

The measurement has been performed using the —)■ 7 r^ 0 (zx’''/i“) as a control 

channel. The only peaking background contributing is that originated from —)■ 
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Figure 4: Observed (solid curve) and expected (dashed curve) CLs values as a function of B[D^ —>• 
7r+7r“^+/r“). The green (yellow) shaded area contains 68.3% and 95.5% of the results of the analysis 
on experiments simulated with no signal. The upper limits at the 90(95)% CL are indicated by the 
dashed (solid) line. 

7 r'''7r’''7r“). Candidates have been selected via a multivariate analysis based on BDT 
(the input variables used are of decay vertex and high distance,p and pt of 

all tracks, IP of all tracks, maximum distance of closest approach between all pairs 
of tracks in the candidate D decay) and some requirements on particle identification 
variables. 

The yields have been obtained with a binned maximum likelihood £t. Details are 
shown in Figj^in m(/i+/r“) bins. 

The results extrapolated to the full dimuon mass spectrum are: 

B{D+ 7r+p+/i-) < 7.3(8.3) ■ IQ-® at 90%(95%) CL 

i3(D+ ^ 7r+/i+/i-) < 4.1(4.8) ■ 10“^ at 90%(95%) CL 

i3(D+ ^ 7r-p+/i+) < 2.2(2.5) ■ 10"® at 90%(95%) CL 

B{Dt 7r-/i+/i+) < 1.2(1.4) ■ 10-^ at 90%(95%) CL 

and represent an improvement by a factor 50 with respect to previous limits, but 
are still one order of magnitude above several NP predictions on —)■ ■ HI. 

6 Future prospects 

The strictest limit on the results currently obtained on charm rare decays field is 
certainly the statistic available. 

As it is well known, LHC is imminently starting with Run II data taking, which 
in the specific case of LHCb will enhance the energy in the centre of mass from 8 to 
13 TeV. At the end of Run II, LHCb will have collected data corresponding to an 
integrated luminosity of 8 fb~^. Assuming that efficiencies and signal-to-background 
ratio will stay unvaried, it is possible to predict an improvement on current limits by 
simply scaling them according to the increased statistics. The improvements due to 
the data collected during Run H should be limited (Tab.l). 
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Figure 5: Invariant mass distributions for: —>■ 7r+/r+/i candidates in the ve ) 

bins: a) low-m(^+/r“) [250-525]MeF/c^, b) rj [525-565]MeF/c^, c) p/w [565-850]MeF/c^, d) (j) 
[850-1250]MeI^/c^, e) high-m(p+p“) [1250-2000]Mey/c^; —)■ 7r“p+p+ candidates in the ve 

m(p+7r-) bins: a) 250-1140 MeV/c^, b) 1140-1340 MeV/c^, c) 1340-1550 MeV/c^, d) 1540-2000 
MeV/c^. 

The data are shown as points (black)and the total PDF (dark blue line) is overlaid. The components 
of the t are also shown: the signal (light green line), the peaking background (solid area) and the 
non-peaking background (dashed line). 
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Figure 6: Observed (solid curve) and expected (dashed curve) CLs values as a function of 
{B)iD+ —>■ TT+p+p ). The green (yellow) shaded area contains the Tier (±2cr) interval of pos¬ 
sible results compatible with the expected value if only background is observed. The upper limits 
at the 90% (95%) CL are indicated by the dashed (solid) line. 
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Mode 

Run I 

Run II 

Upgrade 

—)■ hh' fji~ 

few 10“^ 

fewer 10“^ 

o 

1 

00 

^ /i+/x- 

few 10-9 

fewer lO-^ 

O 

1 

o 

—)■ 7r+/i+/i“ 

few 10-^ 

fewer 10“® 

10-9 

D+ K+fi+fi- 

few 10“^ 

fewer 10“^ 

o 

1 

00 

A —)■ p/r/i 

few 10“^ 

fewer 10“^ 

O 

1 

00 

—)■ e/i 

few 10-^ 

fewer 10“® 

10-9 

(TAcp{D^ 07 ) 

~ 10% 

~5% 



Table 1: Prediction on upper limits on branching fractios,based on Run I results 
(second column). 


Mode 

Run II 

Upgrade 

—)■ TT+zi+p- 

0.6%(30000 ev.) 

0.2%(300000 events) 

—)• 7r~^n~ 

3%(1500 ev.) 

1%(15000 events) 

—)■ R-TT+zi+p- 

1%(10000 ev.) 

0.3%(100000 events) 

K+n-fi+fi- 

40%(30 ev.) 

12%(300 events) 

K+K-^+ijl- 

11%(150 events) 

4% (1500 events) 


Table 2: Prediction on sensitivities to asymmetries in multi-body rare decays. 


After the Upgrade phase, with a/s = 14 TeV, with a total amount of 50 fb~^, one 
more order of magnitude could be reachable (Tab.l), allowing to conhrm or discredit 
some NP models predictions which still lay between current limits and the SM theory. 

As written in the introduction, there is a stimulating held which is still to be 
investigated, and represents a great chance to test NP theories and SM too. Indeed, 
the multibody charm decays allow to explore CP [15] , forward-backward ng-na and 
also T-odd asymmetries [18] helds. Enhancements of 0(1%) asymmetries could be 
measurable in —)■ 7r’''/i/r, while 0(5%) may be detectable with —)■ 

(Tab.2). 

However, these expectations could improve signihcantly under the upgrade condi¬ 
tions due to the contribution of improvements in the analyses, combination of modes, 
which might matter more than individual sensitivities, and an offline reconstrucion 
quality available in a fully software trigger, with an improvement of efficiencies that 
could reach a factor 3. 

7 Conclusion 

The results presented are all best world limits. 
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An update of —)■ is currently in progress and the results obtained in 

the search for the Lepton Flavour Violating —)■ e^/i^ will become public soon. 

Another measurement, the —)■ fi~ branching fraction should very soon 

be public. It will serve as a normalisation mode in the future analyses ssearchinf for 
NP in —)■ decays, that are presently on going. 

In the next decade, the available datasets will be multiplied by two orders of 
magnitude. It will make more sophisticated measurements, like CP or angular asym¬ 
metries, possible. 
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